Accepted for publication in The Astrophysical Journal Letters 

Preprint typeset using L^-T^X style cmulateapj v. 5/2/11 



(N 

o 

Oh: 



O 

u 

6 

in 

o3 



Ha EQUIVALENT WIDTHS FROM THE 3D-HST SURVEY: EVOLUTION WITH REDSHIFT AND 

DEPENDENCE ON STELLAR MASS 

Mattia Fumagalli 1 , Shannon G. Patel 1 , Marijn Franx 1 , Gabriel Brammer 2 , Pieter van Dokkum 3 , Elisabete da 
Cunha 5 , Mariska Kriek 4 , Britt Lundgren 3 , Ivelina Momcheva 3 , Hans-Walter Rix 5 , Rasper B. Schmidt 5 , 
Rosalind E. Skelton 3 , Katherine E. Whitaker 3 , Ivo Labbe 1 , Erica Nelson 3 

Accepted for publication in The Astrophysical Journal Letters 

ABSTRACT 

We investigate the evolution of the Ha equivalent width, EW(Ha), with redshift and its dependence 
on stellar mass, using the first data from the 3D-HST survey, a large spectroscopic Treasury program 
with the HST-WFC3. Combining our Ha measurements of 854 galaxies at 0.8 < z < 1.5 with those 
of ground based surveys at lower and higher redshift, we can consistently determine the evolution of 
the EW(Ha) distribution from z=0 to z=2.2. We find that at all masses the characteristic EW(Ha) 
is decreasing towards the present epoch, and that at each redshift the EW(Ha) is lower for high-mass 
galaxies. We find EW(Ha) ~ (1 + z) 18 with little mass dependence. Qualitatively, this measurement 
is a model- independent confirmation of the evolution of star forming galaxies with redshift. A quan- 
titative conversion of EW(Ha) to sSFR (specific star-formation rate) is model dependent, because of 
differential reddening corrections between the continuum and the Balmer lines. The observed EW(Ha) 
can be reproduced with the characteristic evolutionary history for galaxies, whose star formation rises 
with cosmic time to z ~ 2.5 and then decreases to z = 0. This implies that EW(Ha) rises to 400 A 
at z = 8. The sSFR evolves faster than EW(Ha), as the mass-to-light ratio also evolves with redshift. 
We find that the sSFR evolves as (1 + z) 3 - 2 , nearly independent of mass, consistent with previous 
reddening insensitive estimates. We confirm previous results that the observed slope of the sSFR-z 
relation is steeper than the one predicted by models, but models and observations agree in finding 
little mass dependence. 

Subject headings: galaxies: evolution galaxies: formation galaxies: high-redshift 
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1. INTRODUCTION 

Several studies have combined different star formation 
indicators in order to study the evolution of star-forming 
galaxies (SFGs) with redshift. At a given redshift low 
mass galaxies typically form more stars per unit mass 
(i.e., specific star- formation rate, sSFR) than more mas- 
sive galaxies (Juneau et al. 2005, Zheng et al. 2007, 
Damen et al. 2009). In addition the sSFR of galaxies 
with the same mass increases at higher redshift. How- 
ever, semi-analytical models and observations are at odds 
with regards to the rate of decline of the sSFR towards 
low-redshift (Damen et al. 2009, Guo et al. 2010). 

One of the main observational caveats is that most of 
the studies covering a wide redshift range use diverse 
SFR indicators (such as UV, IR, [Oil], Ha, SED fitting). 
This is a consequence of the fact that it is difficult to use 
the same indicator over a wide range of redshifts. One 
therefore has to rely on various conversion factors, often 
intercalibrated at z — and re-applied at higher redshift. 

A well-calibrated standard indicator of the SFR is the 
Ha luminosity (Kennicutt, 1998). However, Ha is shifted 
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into the infrared at z > 0.5, and it is difficult to mea- 
sure due to the limitations of ground-based near-IR spec- 
troscopy. Comparing measures of Ha at different red- 
shifts has therefore been a challenge. Most of the Ha 
studies at high redshift are based on narrow-band pho- 
tometry (e.g. the HiZELS survey, Geach et al. 2008). 

The 3D-HST survey (Brammer et al., 2012) provides a 
large sample of rest-frame optical spectra with the WFC3 
grism, which includes the Ha emission in the redshift 
range 0.8 < z < 1.5. Taking advantage of the first data 
from the survey (45% of the final survey products) we 
investigate for the first time the star formation history 
(SFH) of the Universe with Ha spectroscopy, using a 
consistent SFR indicator over a wide redshift range. 

We evaluate the dependence of the Ha Equivalent 
Width, EW(Ha), on stellar mass (M*) and redshift (up 
to z ~ 2), comparing the 3D-HST data with other sur- 
veys in mass selected samples with M* > 1O 1O M0. Since 
EW(Ha) is defined as the ratio of the Ha luminosity to 
the underlying stellar continuum, it represents a measure 
of the the current to past average star formation. It is 
therefore a model independent, directly observed proxy 
for sSFR. 

We also derive SFRs from the Ha fluxes. We evaluate 
the mean sSFR in stellar mass bins and study its evo- 
lution with redshift. The slope of the sSFR-z relation 
in different mass bins indicates how fast the star forma- 
tion is quenched in galaxies of various masses. Finally, 
we compare our findings to other studies (both observa- 
tions and models), discussing the physical implications 
and reasons for any disagreements. 
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Fig. 1. — EW(Ha) against mass, for different redshift samples. Vertical lines represent the limiting mass of the analysis. Black symbols 
are objects with Ha detection with S/N > 3 and red arrows represent upper limits. The green diagonal lines represent the detection limit of 
the 3D-HST data. Blue solid lines represent the mean EW(Hq) of detected SFGs, in 0.5 dex mass bins. Red solid lines represent the mean 
EW(Ha) of all galaxies, assuming EW(Ha)=0 for non-detected objects. Errors of the means are computed with a bootstrap approach. At 
each redshift higher mass galaxies have lower EW(Ha) than less massive objects. 



2. DATA 
2.1. 3D-HST 

We select the sample from the first available 3D-HST 
data. They include 25 pointings in the COSMOS field, 6 
in GOODS-South, 12 in AEGIS, 28 in GOODS-Nortlfl 
Spectra have been extracted with the aXe code (Kum- 
mel et al., 2009). Redshifts have been measured via the 
combined photometric and spectroscopic information us- 
ing a modified version of the EAZY code (Brammer, 
van Dokkum, Coppi, 2008), as shown in Brammer et 
al. (2012). Stellar masses were determined using the 
FAST code by Kriek et al. (2009), using Bruzual & 
Chariot (2003) models and assuming a Chabrier (2003) 
IMF. The FAST fitting procedure relies on photometry 
from the NMBS catalogue (Whitaker et al. 2011) for the 
COSMOS and AEGIS fields, the MODS catalog for the 
GOODS-N (Kajisawaet al. 2009) and the FIREWORKS 
catalogue for GOODS-S (Wuyts et al. 2008). 

The mass completeness limit is log(M*/M Q ) > 10 at 
z=1.5 (Wake et al. 2011, Kajisawa et al. 2010). We 
select objects with 0.8 < z < 1.5; this resulted in a 
sample of 2121 galaxies with log(M*/M ) > 10. 

In slitless spectroscopy, spectra can be contaminated 

6 From program GO-11600 (PI: B. Weiner) 



by overlapping spectra of neighboring galaxies. The aXe 
package provides a quantitative estimate of the contam- 
ination as a function of wavelength, which can be sub- 
tracted from the spectra. We conservatively use spectra 
where the average contribution of contaminants is lower 
than 10% and for which more than 75% of the spectrum 
falls on the detector. After this selection we have 854 
objects in the redshift range 0.8 < z < 1.5 (40% of the 
objects). The final sample is not biased with respect to 
the mass relative to the parent sample. 

Line fluxes and EWs were measured as follows Q We 
fit the ID spectra with a gaussian profile plus a linear 
continuum in the region where Ha is expected to lie. 
We subtract the continuum from the fit and measure 
the residual flux within 3er from the line center of the 
gaussian. Errors are evaluated including the contribution 
from the error on the continuum. We distinguish between 
detections and non-detections of Ha with a S/N thresh- 
old of 3. The typical 3tr detection limit corresponds to 
SFR = 2.8M Q yr- 1 at 2=1.5 (Equation 2). 

Due to the low resolution of the WFC3 grism, the Ha 
and [Nil] lines are blended. In this work EW(Ha) there- 
fore includes the contribution from [Nil]. For the other 

7 Through the entire paper the quoted EWs are rest-frame val- 
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Fig. 2. — Evolution of EW(Ha) (top) and sSFR(Ha) (bottom) with redshift, in different mass bins, for SFGs (left) and all objects (right). 
Errors on the average EWs have been evaluating through bootstrapping. Dotted lines are the best fit power laws EW(z) ~ (1 + z) p . At 
fixed mass the average EW(H«) and sSFR(Ha) increase with redshift, with a power law of EW(Hct) ~ (1 + 2) 18 and sSFR(Ha) ~ (l + z) 3,3 
with little mass dependence. The effect of a luminosity dependent dust correction (Garn et al. 2010) correction is shown by the right axis. 
The effect of A(Hc«)=l is shown by the black arrow. 



datasets, which have higher spectral resolution, we com- 
bine Ha and [Nil] for consistency with 3D-HST. 

2.2. SDSS 

We retrieve masses and EW(Hq) for the SDSS galaxies 
from the MPA- JHU catalogue of the SDSS-DR7. Masses 
are computed based on fits to the photometry, follow- 
ing Kauffmann et al. (2003) and Salim et al. (2007). 
At redshift 0.03 < z < 0.06, for masses higher than 
M* = 10 10 Af Q , the SDSS sample is spectroscopically 



complete in stellar mass (Jarle Brinchmann, private com- 
munication). We consider as detections only measure- 
ments greater than 3A, as the ones with EW < 3A are 
affected by uncertainties in the stellar continuum sub- 
traction (Jarle Brinchmann, private communication). 

In the redshift range 0.03 < z < 0.06, the spectroscopic 
fiber of SDSS does not cover the entirety of most galaxies. 
As a consequence sSFR are evaluated with emission line 
fluxes and masses from the fiber alone. 
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2.3. VVDS 

The VIMOS VLT Deep Survey (VVDS, Le Fevre et 
al. 2005) is a wide optical selected survey of distant 
galaxies. Ha is covered by the VIMOS spectrograph at 
0.0 < z < 0.4. 

Lamareille et al. (2009) released a catalog of 20,000 
galaxies with line measurements, complete down to M* = 
10 9 5 M Q at z = 0.5 . Masses are retrieved from the 
VVDS catalog; they have been computed though a 
Bayesian approach based on photometry (equivalent to 
Kauffmann et al. 2003 and Tremonti et al. 2004), and 
are relative to a Chabrier IMF. We select a sample with 
redshift 0.2 < z < 0.4 and M* > 10 10 M Q containing 
741 objects, of which 477 (64 %) have an Ha measure- 
ment with S/N > 3. The percentage of Ha detected 
objects drops to 32% at masses M* > 10 110 M o . 

2.4. High redshift data 

Erb et al. (2006) published EW(Ha) for galaxies se- 
lected with the BX criterion (Steidel et al. 2004), tar- 
geting SFGs at redshift 2.0 < z < 2.5. We evaluate the 
completeness of the sample as follows. From the FIRE- 
WORKS catalogue (Wuyts et al. 2008) we reconstruct 
the BX selection and evaluate the fraction of objects with 
spectroscopically confirmed redshift 2.0 < z < 2.5 that 
fall in the BX selection. Percentages are 44%, 32% and 
27% for mass limited samples with log(M*/M )) = 10.0- 
10.5, log(M*/M ) = 10.5-11, log(M*/M ) > 11.0. 

3. THE EW(Ho) - MASS RELATION 

We first study how EW(Ha) depends on stellar mass 
in each available data set. The 3D-HST sample has 
been divided in two redshift bins, 0.8 < z < 1.1 and 
1.1 < z < 1.5. The results are shown in Figure [1] At 
each redshift, highest mass galaxies have lower EW(Ha). 
Note however, that there is a large scatter in the relation. 
We quantify the trend in the following way: we determine 
the average EW(Ha) in three 0.5 dex wide mass bins 
(10.0 < log(M*/M < 10.5, 10.5 < log(M*/M ) < 11, 
log(M*/M Q ) > 11.0) and evaluated its error through 
bootstrapping the sample. The mean EW(Ha) in a given 
mass bin is obtained in two ways: (1) using only de- 
tected, highly SFGs (blue lines in Figure Q]) and (2) us- 
ing all galaxies but assigning EW(Ha)=0 to the objects 
detected in Ha with S/N < 3 (red lines in Figure [TJ. For 
the z = 2.2 data, we use the the FIREWORKS catalog 
to establish the fraction of galaxies excluded by the BX 
selection and therefore give an estimation of EW(Ha) for 
all galaxies. 

Using either method we find that an EW(Ha)-mass re- 
lation is in place at each redshift, not just for strongly 
star forming objects but also for the entire galaxy pop- 
ulation. Galaxies in the lowest mass bin (10.0 < 
log(M»/M Q ) < 10.5) have on average an EW(Hq) which 
is 5 times higher than galaxies in the highest mass bin 
(log(M»/M ) > 11.0). 

We discuss the evolution of the EW(Ha)-mass relation 
with redshift in the next section of the Letter. However, 
it is immediately evident from Figure[T]that the 3D-HST 
survey targets galaxies with EW(Ha) typical of an inter- 
mediate regime between what is seen at z=0 and what is 
seen at higher redshift. In other words the EW(Ha)-mass 
relation seems to rigidly shift towards higher EW(Ha) at 



higher redshifts. 

A considerable fraction of detected galaxies in 3D-HST 
have EW(Ha) > 30A, while in SDSS similar objects are 
extremely rare: 3.8%, 1.4% and 0.4% for increasing mass 
samples. This study can be seen as an extension of the 
findings of van Dokkum et al. (2011), who reported that 
massive galaxies at z > 1 show a wider range of EW(Ha) 
compared to galaxies in the local Universe. Following 
this trend with redshift, in the z > 2 bin we find typical 
EW(Ha) of 150 A for SFGs. Such high values of EW(Ha) 
represent just 11% of the 3D-HST sample (14.6%, 8.1% 
and 4.7% respectively for increasing mass bins). 

4. THE EVOLUTION OF EW(Ha) WITH REDSHIFT 

The evolution of EW(Ha) with redshift can be seen as 
an observational (i.e. model-independent) proxy for the 
sSFR-z relation. Figure[2](top panels) shows the redshift 
evolution of the average EW(Ha) in different mass bins 
for the detected SFGs (top left panel) and for all galaxies 
(top right panel) . A substantial increase of the EW(Ha) 
is seen at higher redshifts in both samples. We therefore 
infer that evolution of the SFR is not a byproduct of 
selection effects from different SFR indicators. At 0.8 < 
z < 1.5, a galaxy has on average an EW(Ha) that is 3-4 
times higher than that of an object of comparable mass 
in the local universe. For each mass bin we parametrize 
the redshift evolution of the EW(Ha) as follows: 

EW(Ha)(z)~Ax(l + zf (1) 

The coefficient p has an average value of 1.8, with little 
dependence on mass (best fit values are listed in Table 
Q}. As can be deduced from Table 1 there may be a 
weak mass dependence such that the relations steepen 
with mass; however, the difference between the slopes 
at the lowest and highest mass bin is not statistically 
significant. 

This indicates that the decrease of EW(Ha) happens 
at the same rate for all galaxies irrespectively of their 
masses. As seen in the right panel of Figure 2, the ad- 
dition of non-SFGs amounts to a negative vertical shift 
in the EW(Ha) but not to a change in the slope of the 
relation. 

An uncertainty is the effect of dust on the EW(Ha). 
Without more measurements we cannot state what ef- 
fect dust has, and in literature there is disagreement on 
the relative extinction suffered by the nebular emission 
lines and the underlying stellar continuum (Calzetti et 
al. 2000, Erb et al. 2006, Wuyts et al. 2011). However, 
the data motivated model described in Section 6 suggests 
that dust has a mild effect on the EW(Ha). 

5. THE SSFR(Ha) - MASS RELATION AND ITS 
EVOLUTION WITH REDSHIFT 

The EW(Ha) has the advantage that it is a direct ob- 
servable, but it is more difficult to interpret than a more 
physical quantity like the specific star formation rate. 
The latter can only be derived with a proper extinction 
correction for Ha. We lack this information, as we do not 
have a proper Balmer decrement measurement. In the 
following the briefly explore the specific star formation 
rate (sSFR) evolution implied by assuming no extinction 
and later discuss the effect of a dust correction to the 
measured slopes of the sSFR-z relation. 
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Fig. 3. — Mean values of sSFR(Ha) in 0.5 dex mass bins at 
various redshift for SDSS, VVDS, 3D-HST and from Erb et al. 
(2006). At each redshift more massive galaxies have less sSFR(Ha) 
than less massive ones. The effect of a luminosity dependent dust 
correction (Garn et al. 2010) correction is shown by the right axis. 
The effect of A(Hct)=l is shown by the black arrow. 

The SFR is derived from the Ha flu30 using Kcnnicutt 
(1998): 

SFRtHaOlMoyr- 1 ] = 7.9x 1CT 42 xL(Ha)[erg/s] x lCT ' 24 

(2) 

where the 10 24 factor accounts for a conversion to the 
Chabrier IMF, from Salpeter (as in Muzzin et al. 2010). 

Figure[3]shows the mean value of sSFR(Ha) in different 
stellar mass bins, at different redshifts. In each redshift 
bin higher mass galaxies have lower sSFR(Ha). 

In Figure [5] (bottom panels) we show the redshift evo- 
lution of the average sSFR(Ha) in different mass bins 
for detected SFGs (bottom left panel) and for all galax- 
ies (bottom right panel). A typical galaxy at z=1.5 has 
as sSFR(Ha) 15-20 times higher than a galaxy of the 
same mass at z=0. In each mass bin we fit the evolution 
of the sSFR in redshift with a power law: 



sSFR(z) ~ B M x (1 + zf 



(3) 



obtaining a value of n — 3.2 ± 0.1. As can be deduced 
from Table 1 there may be a weak mass dependence such 
that the relations steepen with mass; however, the slopes 
at the lowest and highest mass bin are not statistically 
different. 

The sSFR-z relation is steeper than the EW-z relation, 
because of the additional evolution of the M/L ratio: 

EW/sSFR - L(Ha)/L R x M*/(K*L(Ha)) - M/L R (4) 

where K is the conversion factor in Equation 2 and Lr 
is the R-band luminosity. 

Implementing a luminosity dependent dust correction 
for Ha (Garn et al. 2010, shown with the right axis in 
Figure [2j bottom right) would increase the value of n to 
3.7 ±0.1. However, several studies (Sobral et al. 2012, 

8 We assume a [NII]/(Hq+[NII]) ratio of 0.25 for the 3D-HST 
sources. 



Dominguez et al. 2012, Momcheva et al., submitted) 
have indicated that a better indicator for the Ha extinc- 
tion at different redshifts is the stellar mass, and that 
the Ha extinction depends strongly on mass but little on 
redshift (at constant mass). The Garn & Best 2010 re- 
lation gives median A(Ha) of 1, 1.5 and 1.7 mag for the 
increasing mass bins in this study. A mass-dependent 
dust correction impacts the normalization of the sSFR 
but not the slope n. 

Our implied evolution of the sSFR compares well to 
results from literature. For example, Damen et al. (2009) 
found n = 4 ± 1 based on UV + IR inferred sSFRs, and 
Karim et al. (2011) found n = 3.50 ± 0.02 for SFGs and 
n — 4.29 ± 0.03 for all galaxies, based on stacked radio 
imaging. 

All results indicate an evolution which is steeper with 
redshift than semianalytical models (Guo & White 2008, 
Guo et al. 2011), who find slopes close to n=2.5. All 
studies find that the slope does not depend on the stellar 
mass out so z=2. In short, our results are consistent with 
previous determinations. 

6. LINKING THE CHARACTERISTIC SFH OF GALAXIES 
AND EW(Ha) 

We compare the observed evolution of EW(Ha) to 
what might be expected from other observations. We 
construct the typical SFH of a galaxy with mass ~ 
10 11 M Q at z=0. As a starting point, we assume that 
the cumulative number density remains constant with 
redshift (similar to van Dokkum et al. 2010, Papovich et 
al. 2011, Patel et al. 2012). We use the mass functions 
of Marchesini et al. (2009) and Papovich et al. (2011), 
and we show the resulting mass evolution in Figure Hp. 
We determine the SFR at these masses from Damen et 
al. (2009), Papovich et al. (2011) and Smit et al. (2012), 
and we fit a simple curve to these values (indicated by 
the curve in Figure 0^,) . This evolutionary history repro- 
duces the mass evolution well (Figure 0J:) . 

Next we calculate the implied EW(Ha): L(Ha) is de- 
rived using Equation 2, and Bruzual & Chariot 2003 
models are used to calculate the stellar continuum, as- 
suming solar metallicity (Figure HJd). The predicted 
EW(Ha) rises monotonically to high redshift, reaching 
400A at z=8 (Figure Hi). The predicted EW(Ha) corre- 
sponds surprisingly well to the observed EW(Ha) Even 
the z=4 detections by Shim et al. (2011) based on broad- 
band IRAC photometry are consistent within the errors. 
Apparently, our simple method produces a robust pre- 
diction of the evolution of EW(Ha). We note that the 
implied sSFR (Figure HJj) is higher than expected from 
straight measurement of L(Ha), consistently with signif- 
icant dust extinction. One magnitude of extinction for 
Ha is needed to reconcile this discrepancy. 

It is remarkable that our prediction worked well for 
EW(Ha): the average evolution of galaxies was derived 
from the evolution of the mass function and SFR, which 
carry significant (systematic) uncertainties when derived 
from observations; whereas the EW(Ha) is a direct ob- 
servable. 

7. CONCLUSIONS 

We have used the 3D-HST survey to measure the evo- 
lution of the EW(Ha) from z=0 to z=2. We show that 
the EW(Ha) evolves strongly with redshift, at a con- 
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Fig. 4. — Comparison of observed EW(Ho) with predictions from a simple observational supported model, at different redshifts. (a) 
Input SFH, and Ha Luminosity (b) Luminosities at 6563 A, from the Bruzual & Chariot 2003 code, (c) Mass growth, (d) Evolution of 
EW(Ha) with redshift. (e) Evolution of sSFR with redshift. Data points are mean EW(Ha)/sSFR(Ha) of observed galaxies with mass in 
a 0.3 dex bin around the typical mass of the model at a given redshift. 



stant mass, like (1 + z) 1 ' 8 . The evolution is independent 
of stellar mass. The equivalent width goes down with 
mass (at constant redshift). The increase with redshift 
demonstrates the strong evolution of star forming galax- 
ies, using a consistent and completely model independent 
indicator. We explore briefly the implied sSFR evolution, 
ignoring dust extinction. We find that the evolution with 
redshift is strong (sSFR ~ (1 + z) 3 2 ). This stronger evo- 
lution is expected as the mass-to-light ratio of galaxies 
evolves with time, and this enters the correction from 
EW to sSFR. The increase with redshift is faster that 
predicted by semi-analytical models (e.g., Guo & White 
2008), consistent with earlier results. 

We construct the characteristic SFH of a 10 11 M Q 
galaxy. This simple history reproduces the observed evo- 
lution of the EW(Ha) to z=2.5, and even to z—A. It 
implies that the EW(Ha) continue to increase to higher 
redshifts, up to 400 A at z=8. This has a significant im- 



pact for the photometry and spectroscopy of these high 
redshift sources. 

The study can be expanded in the future when the 
entire 3D-HST survey will be available, doubling the 
sample and including the ACS grism. In addition to in- 
creased statistics, the ACS grism will allow evaluation of 
the Balmer decrement and therefore a precise dust cor- 
rected evaluation of SFR. Moreover, a statistically signif- 
icant Ha sample at z ~ 1 will be central to understand 
the composition, the scatter and the physical origin of 
the so called 'star- forming- main sequence'. 
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TABLE 1 

Slopes of EW-z and sSFR-z Relation 



log(M»/M Q ) EW(det) EW(all) sSFR(dct) sSFR(all) 



10.0-10.5 1.79 ±0.18 1.52 ±0.21 3.32 ± 0.08 3.06 ±0.13 

10.5-11.0 1.89 ±0.20 1.75 ±0.13 3.18 ±0.09 3.11 ±0.18 

11.0-11.5 2.21 ±0.22 2.12 ±0.43 3.50 ±0.12 3.45 ± 0.26 



